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Abstract 

A Thomas-Fermi model of a spherical shell of positive charge is 
investigated, under various boundary conditions. The electron distri- 
bution and the ionization charge are given particular attention. 
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The advent of the fullerene molecule |1|] and the metallic clusters|2[] have 
called the attention upon the stability of new atomic micro-objects. One of 
them is described in this paper. As it is known, the fullerene molecule Cqq 
consists of 60 carbon atoms (C), arranged (in pentagons and hexagons) on the 
surface of a sphere of radius ~ 3.5 A. Certain clusters, made of (relatively) a 
small number of atoms (as, for example, alkali atoms), have been identified in 
solid-state matrices, particularly in the octahedral interstices of face-centered 
cubic fullerites. They acquire regular geometrical shapes, like tetrahedrons, 
or cubes, the latter being sometimes centered. We investigate herein the 
stability of a Thomas-Fermi model suggested by these micro-objects. 

As it is known, the Thomas-Fermi model starts with free electrons and 
assumes that, due to the Pauli exclusion principle, there is a certain scale 
length over which their number, and their effective interaction, vary slowly. 
It is a quasi-classical theory, and its validity is ensured by the number of elec- 
trons being much greater than unity. Of course, at infinitely long distances 
the theory is not valid, nor at very short distances where the Coulomb po- 
tential is singular. In the case of an atom with the atomic number Z we 
know that the Thomas-Fermi model is not valid for distances shorter than 
~ 1/Z. As in our case of the micro-objects presented above we are interested 
in distances of the order of the size of large molecules, or of the order of the 
solid-state distances, we may also neglect the variations over atomic scale 
lengths, i.e. over the Bohr radius an = h 2 /me 2 = 0.53 A (where m is the 
electron mass, e is the electron charge and h is the Planck's constant). This 
allows one to treat the atoms in the above micro-objects as being uniformly 
distributed over a spherical shell of radius R. In addition, one may assume 
that the centre of the sphere is occupied by a nucleus of positive charge z Q . 

A positive charge z uniformly distributed over the surface of a sphere of 
radius R, and a central positive charge zo, create an electric potential 

V x (r) = z/R+ z /r , r<R , (1) 

and 

V 2 (r) = (z + z )/r , r> R . (2) 

An electric field 

E s = (z/2 + z )/R 2 (3) 
acts outwardly on the spherical surface, which tends to blow the sphere up; 
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and the electrostatic energy of the object is 



So = z (z/2 + z ) /R . (4) 

The density of free electrons n is related by the Fermi wavevector k F 
through n = k F /3ir 2 ; the local energy of the electrons should be a constant, 
for equilibrium, 

h 2 

- — kp — etp — eV = const , (5) 

where cp is the electrostatic potential of the electron distribution, and V = 
Vi j2 for r < R and r > R, respectively. We shall assume that the electron 
distribution extends to infinity, in which case const = in (5). We shall also 
use the atomic units a# and e 2 /a# = 27.2 eV, which together with h — 1 
give m = 1 and e 2 = 1. Then, from (5), we have 

n=^ + V)i (6) 
and the Poisson equation reads 

+ =47m = ^( V + V)l (7) 

for r ^ R; remark that AV = 0. Introducing the reduced variable x = r/R 
and 

97T 2 Y /oN 
V + V = j— , 8 

^ 128i? 4 x ' V ; 

we get from (7) the Thomas-Fermi equation 

In 3 

=X 2 • (9) 



Since n is a continuous function ip and its two first derivatives are continuous; 
therefore x has a slope discontinuity (but itself is a continuous function), 
exactly as the derivative of V. Defining xi — X f° r x < 1 and Xi = X f° r 
x > 1, we have therefore 

Xi(l)=X 2 (l) • (10) 

The number jV(x) of electrons inside a sphere of radius x is easily obtained 
from 

N(x) = AttR* J o dx x 2 n = ^ dx (x X - x) , (H) 
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whence 



and 



N2{x) = -^{xx2-X2+[xi{l)-X2^)}+xM} , x>l . (13) 



9tt 2 



128i? 3 



xx'i ~ Xi + Xi(0) 



x < 1 



(12) 



128i? 3 

We remark that N(x) is continuous at x — 1, and 

9tt 2 



iV(l) 



128i? 3 



xi(i)-xi(i)+xi(o) 



(14) 



As the system extends to infinity we have to assume that X2(oo) = (together 
with its derivatives), so that the total number of electrons is given by 



N = 



9tt 2 



128.R 3 



xi(i)-x' 2 (i)+xi(o) 



(15) 



The electric field Ei(x) at x < 1, and the total charge qi(x) inside the 
sphere of radius x < 1, are obtained easily from the Gauss' law, 



d 



9tt 2 



4vrr 2 J B 1 = -AttRx 2 ^- fa + V l ) = ^-^-^ (xi - x X [) = ^qi 
whence 



Qi 



and 



9^ ( _ >\ 
128i? 3 V Xl XX V ' 

Qi(x) 



E,{x) = 



R 2 x 2 



(16) 

(17) 
(18) 



Since q± = z — N±, we obtain from (12) and (17) 

9tt 2 



rXi(O) • 



128i? 3 

Similarly, the electric field at x > 1 is E 2 (x) = q 2 (x)/R 2 x 2 and the charge 

9tt 2 



(19) 



12 



128R 3 



(X2 - x X ' 2 ) 



(20) 



4 



whence, by using (13) and (19), we get 

^=l^(xi(l)-xi(l)) • (21) 

We remark that (21) expresses exactly the jump in the slope of x and of V at 
x = 1, as we said above; and the total charge q = ^(oo) = 0,1 e. the infinite 
system is neutral. 

The electrons act on the shell with an electric field E e i given by 



F , - I „- 9?r 2 9 (xA I ,m I „- 

^el dr \r=R 128 R5 dx\x) \ x=1 ^ dr \ r=R ~ 

fx'i(i)-xi(i)l 



9tt^ 
128i? 5 



R 2 



(22) 



or, equivalently, 

E e i 



X 2 (l) - X 2 (1)J - , (23) 



128i? 5 

if we use (10) and (21). From (14) and (19) we find that 

E el = -N(l)/R 2 , (24) 

as expected. In order to have the equilibrium of the shell this field must 
compensate the field E s given by (3), i.e. 

N(l)=z/2 + z ; (25) 

or, using (14), (19) and (21), 

2xi(l)=xi(l)+X 2 (l) • (26) 

Equations (19) and (21) may be viewed as giving the parameters zqR 3 and 
zR 3 , respectively; therefore, we have to solve the Thomas-Fermi equation 
(9) under the rather natural conditions Xi (1) — X2 (1), X2 (00) = and (26). 
There is no such a solution. We have always, in fact, 2%i(l) > x'i(l) + X2(1)> 
which means that N(l) < z/2+zq, i.e. the electrons inside the sphere are not 
numerous enough to ensure the equilibrium; due to their fermionic nature 
they prefer to go outside the sphere where their kinetic energy is lower. 
The infinite Thomas- Fermi "hollow" atom is too "rarefied" to be stable. 
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Obviously, the only way to attain the equilibrium of such a "hollow" atom 
is to embed it into a cage, as, in fact, it is a more realistic case. 

Suppose that we have a spherical cage of radius R > R, where the 
"hollow" atom is introduced. The Thomas- Fermi relationship (5) now reads 



Ti 1 2 

— hp — eip — eV — U = -ey? 
2m 



(27) 



where U is the potential well of the cage and ipo is the chemical potential, 
which must equal the external potential for preventing the flux of electrons 
from either going out or in the cage. Introducing the reduced variable x = 
v/Rq and defining 



up + V + U - v?o 



9tt 2 



X 



128Rk x 



(28) 



(in atomic units) we arrive at the Thomas-Fermi equation (9) with th conti- 
nuity condition (10) at a = R/R . The number of electrons (12) and (13) is 
now given by 



9tt s 



128i?g 



xx'i -Xi + Xi(0) 



x < a 



(29) 



and 

N 2 (x) 



9tt 2 



128#g 



{zx 2 -X2 + a[xi(l)-X 2 (l)] +Xi(0)} , x>a . (30) 



Similar relations (17) — (20) hold now, with R replaced by Ro; while the 
discontinuity condition (21) becomes now 



9tt 2 



128i?g 



(31) 



A similar condition (25) for equilibrium is also obtained. Summing up all 

1 // 3 

these relationships we have to solve the Thomas-Fermi equation x^\ — X 5 
under the following conditions: 



9tt 2 



zo = 



Xi(0) 



128R 3 
Xi(a) = X2(a) 



(32) 
(33) 
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97T 2 

Xi(a)-X2(a) , (34) 



128i2jj 



2xi(a) = a [x'i{a) + x' 2 ( a )J 5 (35) 
in which case the total charge in the cage is given by 

9tt 2 



128.Rjj 



X2(l)-X 2 (l) • (36) 



First, we remark that, as we said above, leaving aside (32) and (34), as 
giving the parameters zoR^ and zRq, requiring X2(l) = and letting R Q go to 
infinity, we have the infinite "hollow" atom discussed previously; and, then, 
it is easily to see that the equilibrium condition (35) is not satisfied, as it 
would require Xi{ a ) = 0? *- e - a vanishing solution. Secondly, we see that if 
we put x 2 (l) — and a = 1 we get again the previous case of an infinite 
Thomas-Fermi atom, which we know that it is unstable; it follows that even 
more unstable will be the "hollow" system with X2(l) = and a < 1, i.e the 
"positive ion". But we remark that this is only a particular cas of a positive 
ion. 

In the remaining of this paper we shall discuss a few types of solutions for 
the Thomas- Fermi equation (9) under the boundary conditions (32) — (35), 
being especially interested in the total charge (36). 

For z = 0, R = 2 A, R = 1.73 A {a = 0.86) and z = 44, as for a 
(tetrahedral) cluster of four sodium atoms, the function x is plotted in Fig.l 
vs x; it corresponds to a total charge q = +2.7. For z = 0, Rq = 3.2 A, 
R = 2.78 A (a = 0.87) and z = 88, as for a (cubic) cluster of eight sodium 
atoms, the function x is shown in Fig.2; the total charge in this case is 
q = -0.1. For z = 11, R = 3.15 A, R = 2.75 A (a = 0.87) and z = 88, as 
for a centered (cubic) cluster of nine sodium atoms, the function x is given in 
Fig.3, for a total charge q = 0.7. We remark that, indeed, the function x has 
a smoother variation with increasing the number of electrons, except for a 
range ~ an around the positions of atomic charges. In addition, we remark 
that the charge q is very sensitive to the input parameters, for increasing 
both the number of electrons and the positive charges. 

We may estimate the energy of the system as follows. The density of 
kinetic electron energy is 

km ~(2n) 3 Jo dk 2 k ~ 10/C ' V 7 > 
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whence the total kinetic energy of electrons 



3 (3-/8) 4 ,1 ^ X 



1 V 5 / 2 



^=^-y dx —> ■ (38) 

The density of potential energy of the electrons is given by 

, rA 3(3tt/32) 3 X 5/2 9tt X 3/2 
^ = -(y, + V + f/)n = , (39) 

whence their total energy 

(3tt/8) 4 r 1 y 5/2 
£ el = E km + E pot = "\^f J Q dx^- N Vo , (40) 

where N = — (q — z — zq) is the total number of electrons. The energy of the 
shell is given by 

£ s = £ -(z + z )U . (41) 

On the other hand we may express the energy of interaction of the shell with 
the electrons in two distinct ways 

Si = z<p(R) = -J dr nV , (42) 
where V is the potential of the shell; from (42) we get 

U -n = ^X2M , (43) 
which together with (40) and (41) allows one to write the total energy as 

£ = £ + £i + £2- [2(z + z ) -q]<p , (44) 

where 

^-^!>^ 

and 

S2 = ~^m {Z + Z ° )X ' 2{1) ■ (46) 



8 



In the first case shown in Fig.l Si = —28, £2 = —42, (x' 2 (1) = 282) while 
the self-energy (4) of the shell is So = 297; for stability, i.e. negative total 
energy, one needs (p > (297 — 28 — 42) /88 ~ 2.6. One can see that, indeed, 
the object is squeezed into the atomic environment, ipo being a measure of 
the variation of the atomic "pseudo-potential" felt by an outside electron 
on its attempts of penetrating the electron cloud of the atomic surrounding; 
these "pseudo-potentials" are potential barriers which confine the clusters. A 
huge "pressure" is exerted by the cluster on its surrounding, which results in 
the deformation of the electronic clouds of the cage walls. Similar values are 
obtained in the other cases, for example, y?o > (742 — 43 — 92) /176 ~ 3.4 
for eight atoms, and y?o > (931 — 155 — 110) /200 ~ 2.4 for nine atoms. 

It might be of interest the variation of q with z. For example, we define 
Z = 128RIz/9tc 2 and Q = l28Rlq/9ir 2 , and solve the equation for various 
values of a. Such a dependence of Q on Z is shown in Fig A for z = and 
a = 0.8, and in Fig.5 for z = and a = 0.9; while in Fig.6 it is shown an 
almost neutral cluster for zq = and a = 0.87. Similar results can also be 
obtained for zo 7^ 0. In Fig.l the variation of q with R Q is shown for zo = 0, 
z = 44 and R — 1.73 A, while in Fig. 8 a similar dependence is included for 
z = 11,2 = 88 and R = 2.75 A. 

In the limit of large number of atoms in the cluster the Thomas-Fermi 
theory is valid. For a finite number n of atoms disposed on a spherical surface 
one may estimate the error in the total charge as follows. From the Poisson 
equation we have that the charge q is proportional to the radial 5ip r and 
angular 5ip a variations of the potential derivatives in the following way: 

q~8<p' r AS + 28<p' a AS , (47) 

where AS is the element of area. Assuming the same variation per unit 
length we find 

Sl ?'a = S<Pr\f^ ■ (48) 

On the other hand, if one neglects the angular variations we have 

g ~(^r) A5 • ( 49 ) 

As these small variations are proportional to the small variations of the dis- 
tance we should also have 

(^;) 2 +2(^:) 2 = (i+87r/n)(^;) 2 = (^;)o , m 
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whence 

J I + 87r/n 

Ag/g~l-^ =4= . (51) 

1 + 2^4n/n 

In our cases of n ~ 4, 8, 9 this error in Ag is about 40%. 
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Figure captions 
Fig.l. 

Function x vs x for z = 0, R = 2 A, R = 1.73 A (a = 0.86) and z = 44, 
and 

a total charge q = +2.7. 
Fig.l. 

Function x vs x for z = 0, i? = 3.2 A, i? = 2.78 A (a = 0.87) and 

2 = 88, 

and a total charge q = —0.1. 
F^.3. 

Function x vs x f° r = H, Ro = 3.15 A, i? = 2.75 A (a = 0.87) and 
^ = 88, 

and a total charge q = 0.7. 
Fi#.4. 

The reduced charge Q vs Z for = and a = 0.8. 
Fig.5. 

The reduced charge Q vs Z for z = and a = 0.9. 
Fig.6. 

The reduced charge Q vs Z for z = and a = 0.87 for an almost neutral 

cluster. 

Fig.l. 

The total charge q vs R for z = 0, z = 44, and i? = 1.73 A. 
Fig.B. 

The total charge q vs R for z — 11, z — 88, and i? = 2.75 A. 
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